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Calculations of the thermal history of the moon were carried out by solving the thermal 
conductivity equation for the case in which the heat sources are the long-lived radioactive 
elements Th, U, and K“. The concentrations of these elements were adjusted to give 
four variations of the heat flow: 1.35 X 10‘* and 0.91 x lO"* cal cm'* s*' (I and I**), 0.61 
X 10'" (variant II, the terrestrial mixture of Lyubimova and Starkova), and 0.236 x lO"” 
(variant III, the chondrite model of Urey and MacDonald). In the same calculations, we 
considered layering of the differentiated material with transport to the surface of the 
radioactive elements after the temperature of the layer rose to 200 K above the melting 
temperature, which is considered in five variants that differ in the amount of transported 
radioactive elements: 100 percent (n = 1), 80 percent n = 0.8), 40 percent, 20 percent 
(n = 0.6, 0.4, 0.2). During fusion the heat capacity and heat conductivity were changed. 

We considered two variants of an initially cold (273 K) and hot (900 K) Moon. Calcula- 
tions show that the interior of the Moon was heated to melting during the first 0.7 to 
2.3 X 10“ years. The maximum fusion involved practically the entire Moon to a distance 
from 15 to 45 km beneath the surface, and started 3.5 to 4.0 x 10’ years ago (I, I**), or 
2.5 to 3.0 X 10’ years ago (II, III) and continued for 1 to 2 X 10* years. Today the Moon is 
cooling. The current thickness of the solid crust is from 150 to 200 km and the heat flow 
exceeds the stationary value 1.5 fold. Apparently the most realistic variant is II (terrestrial 
mixture) for an initially hot Moon, and gives, regardless of the moderate concentration of 
radioactive elements, a heat flow of 0.9 to 0.95 X lO"* cal cm s*‘, which agrees with the 
radioastronomical measurements of Troitsky and Krotikov and with the Apollo 15 data. 


Papers by Urey, MacDonald, Lubimova, 
Levin, Mayeva, Fricker, Reynolds and Sum- 
mers, Iriyama and Shimazu, Anderson and 
Phinney, Hanks and Anderson, MacConnel, 
Toksoz and Solomon, et al. (refs. 1-21) in- 
vestigated the Moon’s thermal history by 
solving the heat conductivity equation with 
given values for the parameters of lunar ma- 
terial and for the initial and boundary con- 
ditions, and with radioactive elements with 
the given concentrations as the heat sources. 
Our papers (refs. 21 and 22) are devoted to 
the same problem. Up to 1966-1967, homo- 
geneous and layered lunar models were con- 


sidered separately. In 1966-1967 we took 
into account in the calculations the differen- 
tiation of lunar material and convection of 
radioactive elements to the surface on melt- 
ing. Working independently, P. Fricker, R. 
Reynolds, and A. Sommers conducted similar 
studies. The consideration of continuous dif- 
ferentiation in the calculation permits one to 
investigate the Moon’s thermal history in 
more detail. In 1971 differentiation was also 
taken into account by Mayeva (ref. 13) and 
in 1973 by Toksoz and Solomon (ref. 21). 

Thus, at the present time computer calcu- 
lation schemes of lunar thermal history are 
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sufficiently fully developed and many new 
results have been obtained. However, inter- 
est in this problem has not yet decreased, but 
rather has increased. This is associated with 
the large uncertainty in the lunar parameters 
that enter into the heat conductivity equa- 
tion. This applies also to the uncertainty in 
the values for density, heat capacity, and 
heat conductivity, and their dependence on 
the state of lunar material (for example, 
during melting) . Also involved is uncertainty 
in the melting curve and the initial condi- 
tions (in particular, the initial temperature) 
and in the Moon’s age. However, the deter- 
mining factor in the Moon’s thermal history 
is the accepted concentration of radioactive 
elements in the Moon and their redistribution 
with time in the lunar interior. 

In the above papers, except for references 
22 and 23, the concentration of radioactive 
elements is taken to be that of chondritic 
meteorites (for example, references 1-5, or 
of the Earth (refs. 7-13), or of a “terrestrial 
mixture,’’ that lead to values of the heat flux 
through the lunar surface on the order of 
(0.3 to 0.6) X 10'® cal cm~^ s"^. At the same 
time radio astronomical investigations by 
V. S. Troitsky and V. D. Krotikov (refs. 24 
and 25) give a flux value equal to 0.85 to 
0.95 X 10“* cal cm“® s"^ which was long 
considered to be rather doubtful but was 
confirmed by the recent analyses of the lunar 
regolith by Apollo 15 (ref. 26). As a result, 
it is necessary to calculate the Moon’s ther- 
mal history for concentrations of the radio- 
active elements that may account for such a 
large heat flux. 

The first attempt at this calculation was 
made in our papers in 1966-1967 (refs. 22 
and 23) where we used four variants of the 
concentration of radioactive elements, two 
of which correspond to concentrations used 
by Levin and Mayeva (ref. 8) (variant “C’’ 
uses the mean terrestrial value of Urey and 
MacDonald). The other two use relatively 
large concentrations which give heat flux 
values close to the radio astronomical val- 
ues. An exponential differentiation of radio- 
active elements was assumed to occur on 
melting. Since the large concentration of ra- 


dioactive elements was close to the surface 
of the Moon (in 10-20 km they weakly in- 
fluence the temperature of the main lunar 
mass, which quickly cooled). As a result, 
the thickness of the solid crust in these 
variants was equal to 600-700 km, i.e., larger 
than for the mean (Levin) and the low 
(Urey, MacDonald) concentrations where 
the thickness was 250 km and 400 km. 

In our studies (refs. 21 and 22), we used 
the exponential model for transport of the 
radioactive elements to the surface on differ- 
entiation which corresponds in some degree 
to the so-called “sudden differentiation.’’ 
This model permits one to determine the 
main physical peculiarities of the thermal 
history related to differentiation ; yet it 
limits the number of variable parameters. In 
this respect, the layered differentiation pro- 
posed by Fricker et al. (ref. 14) appears to 
be closer to the real processes. Therefore, in 
the present paper we have used a layered 
differentiation of radioactive elements on 
melting with different portions of the ele- 
ments (from 20 percent up to 100 percent) 
transported from layer to layer. The concen- 
tration variants of radioactive elements were 
the same as in our 1966-1967 papers (refs. 
22 and 23). 

It should be noted that we are not special- 
ists in geophysics and have no right to dis- 
cuss questions of the Moon’s origin and its 
chemical composition. For this reason, for 
calculations we have taken parameters either 
from the expected values of the heat flux (at 
given concentrations of radioactive elements) 
or from papers by other authors. Thus, the 
initial temperatures (273 K and 900 K) 
were taken on the basis of references 14 and 
16, and the heat capacity on the basis of the 
papers by Mayeva (refs. 11 and 12), where 
a sharp change in the heat capacity at the 
time of melting is considered. Remarks by 
Lubimova (ref. 6) are taken into account 
on the possibility of the flux increasing at 
moderate concentrations of radioactive ele- 
ments by means of the contribution from the 
exciton component of the heat conductivity 
and from amplification of the efficiency of the 
radiational component. 


CALCULATIONS OF THE MOON’S THERMAL HISTORY 


349 


C[r,r(D] 



0.25 - 
4L 

(Ar)2 

4L 

{^T)^ 


o ^ - , at T and T >r,n.„ 

- T, <r 

(T + ^at, 


1. +A2’, 

. +a7’/2, 
<T 


+Ar 


( 1 ) 


f(T) = A/T; A = 10.3 X 10" cinfear '’ ^ ^ >7’in.m.+Ar; 

f(r) = 0.04 X 10" — = const, at <T < (r,„.„.+A7’). 
cm y©&x* 


Consideration is given to the considerable 
idealization, primarily of the layering pro- 
cess, and therefore, the results obtained are 
correct only within the limits of our mathe- 
matical model. 


Statement of the Problem 

As is the custom, the Moon is assumed 
to be a sphere with the radius r = 1735 km, 
of homogeneous density p = 3.34 g/cm®, 
with heat capacity C and the heat conductiv- 
ity K, defined according to Lubimova (ref. 
6) (heat conductivity) and Mayeva (refs. 11 
and 12) (heat capacity) as follows. 

The heat capacity is shown in equation (1). 
L = 100 cal/g is the heat of fusion ; aT = 
200 K; 

r = in cm, T = in K. 

T.nit. melt = 1373 + 500 (1 - rVr.®.) is the 
beginning of melting of the silicate material 
(ref. 3). 

Expression (1) considers the gradual ab- 
sorption of the heat of fusion in the melting 
interval aT (the effective heat capacity in the 
interval AT first increases linearly and then 
decreases) . 

The heat conductivity is given by the expres- 
sion: 


K = f(T) + 


16n®87^ 


In (2) the lattice (phonon), radiational 
(photon) and exciton components of the heat 
conductivity (ref. 6) are taken into account. 

Phonon heat conductivity outside the melt- 
ing interval varies in inverse proportion to 
the temperature; during melting it sharply 
increases (see ref. 27) and remains con- 
stant as shown in the equation above. 
Photon and exciton components of the heat 
conductivity are defined by the well-known 
relations where n® = 3 is the refractive in- 
dex, S = 4.27 X 10 ® cal/cm® grad^ year is 
the Stephan-Boltzmann constant, e = 25 cm'* 
is the absorption coefficient of the material, 
and E = 5.17 X lO-®® cal is the activation * 
energy of excitons. 

For the assumption of a homogeneous dis- 
tribution of heat sources and location of the 
initial coordinate in the center of the sphere, 
the heat conductivity equation has the form: 


pC 


dt r® dr \ dr ) 


+ H, 


(3) 


where H is the heat due to the long-lived 
radioactive elements uranium, thorium, and 
potassium 40. 


H = e ^‘^F{r,t) pai, (4) 

i=i 

where the function F {r,t) defines the distri- 
bution of radioactive elements as a function 
of the radius at different times. Indexes j = 
1,2,3 refer to uranium, thorium, and potas- 


1.76 T® 


rf— 

l\kT 


+ 2 


)V 2 ] 


’The values are taken for pyrolite (a mixture of 
(2) basalt and periodotite) lunar composition, which is 
similar to the Earth’s mantle (see ref. 28). 
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Table 1. — The U, Th, and, K Concentrations in Several Variants of the Bulk 
Lunar Composition and the Stationary Heat Flow Values for These 
Variants 



Concentration of Radioactive 
Elements [g/g] 

Qitat 




Uranium 

Thorium 

Potassium'**" 

cal cm'* S'* 




H ' 1 ° 

H’? -4H*« 

H V 


Basalt 

14% ' 

1 





Peridotite 

Iron 

73% 1 
13% 1 

I 

11.56 X 10-* 

46.7 X 10-* 

18.5 X 10-* 

1.35 X 10- 

(K"/U) = 

1.6 > 

1 





Basalt 
(K“/U) = 

10% j 
1.13 I 

1 

1 

8.60 X 10-’ 

34.0 X 10-* 

9.6 X 10-' 

0.91 X 10- 

Basalt 
Variant by 
Levin 

6% j 
( 

1 

5.20 X 10-* 

21.0 X 10-* 

8.4 X 10-' 

0.61 X 10- 

K"/U = 1.6 J 






Chondrite Model 1 






Variants by 
Urey and 
MacDonald 

1 

) 

i III 

1 

1.16 X 10-* 

4.67 X 10-* 

9.35 X 10-* 

0.236 X 10- 


slum 40, respectively, Xi = 1.54 X 10~^® 
year~S A2 = 5 X lO-^^ year-\ A3 = 5.7 X 
year~^, the heat released by one gram 
of the radioactive element is <»< = 0.805 cal 
g-i year-1, = 0.193, «s = 0.224 and values 
of Hf>, the concentration at the beginning are 
defined on the basis of the expected present 
concentration of H'p according to the for- 
mula 

^ o _ ^ to g ' *°where U = 4.7 X 10® years 

is the age of the Moon. Let us consider four 
variants of the concentrations of radioactive 
elements given in table 1. The concentrations 
of radioactive elements were chosen on the 
basis of the stationary heat fiuxes equal to 
1.35 and 0.91 cal cm-® s-i for variants I 
and and correspond to the radio as- 
tronomical data by Krotikov and Troitsky 
(ref. 24) and to the Apollo 15 data (ref. 26) 
Stationary heat fluxes 0.61 and 0.236 cal 
cm-® s-i for variants II and III coincide 


with one of the variants of Levin and Mayeva 
(ref. 8) (II) and are close to concentrations 
taken by Urey and MacDonald (III). The 
ratio of thorium to uranium is taken equal 
to 4 (see ref. 8) . The uranium concentration 
used in variants I, I^, and II is based on the 
pyrolite lunar composition and the chondrite 
composition in variant III. The assumption 
of the pyrolite lunar composition is based 
on the similarity of the composition of the 
Moon to the Earth’s mantle noted by Ander- 
son (ref. 29) (for example, the mean atomic 
weight of the Earth’s mantle is M = 22.4, 
which refers to eruptive rocks and carbo- 
naceous, ordinary, and enstatite chondrites 
and has the same mean atomic weight as 
lunar material — 22.0). 

The assumption of pyrolite lunar compo- 
sition also is based on the data from Sur- 
veyor V and on the papers by McCrea (ref. 
30) and Reynolds and Summers (ref. 15). 
From all these data it follows that the Moon 
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Rocks : 

Granite 

Basalt 

Eclogite 

Peridotite 

Dunite 

Chondrites 

U X 10-“ 

400 

80 

4.3 

0.6 

0.1 

1.1 

K‘VU 

1.05 

1.13 

1.44 

1.2 

1.2 

9 


consists of silicates (for example, basalt) 
similar to those in the Earth’s mantle.^ The 
above table gives uranium concentrations 
and the ratio of to uranium for different 
rocks. 

From this table and the given stationary 
fluxes we chose the variants I, F"- and II with 
the content of basalt equal to 14 percent, 
10 percent, and 5 percent, respectively; with 
periodotite equal to 73 percent, 77 percent, 
and 82 percent, respectively; and iron, 13 
percent (ref. 15) with an impoverished con- 
tent of radioactive elements (iron has U = 
0.1 X 10-8 g/g, Th = 0.4 X 10-8, and K« = 
0) and the ratio of K*®/U = 1.6, 1.13, and 
1 6, respectively (see refs. 9 and 6). Variant 
III is the chondrite model (Urey and Mac- 
Donald) . Note that variant II corresponds to 
variant “C” with the mean content of radio- 
active elements that of the Earth (see refs. 
31 and 6U. 

As initial conditions we take two variants 
of the “cold” and the “hot” models for the 
formation of the Moon with initial tempera- 
tures of 273° K and 900° K. These are the 
same initial conditions taken by Urey, Mac- 
Donald, Iriyama, Shimazu, and Fricker ^ 
based on the different hypothesis of the 
Moon’s origin and the time of its accumula- 
tion and initial heating, on which we will 
not dwell. 

The boundary condition (r = r„) is de- 
fined by the following relation 



* It does not necessarily follow that the hypothesis 
is true for formation of the Moon by fusion from 
the Earth, since even interstellar material contains 
grains of refractory silicates (see reference 30). 

’ Running ahead, we note that this variant espe- 
cially gives the flux corresponding to Apollo 15 data 
and the radio astronomical data, i.e., evidently, the 
terrestrial composition of the Moon is rather 
probable. 

’ Mayeva (ref. 13) used a parabolic temperature 
distribution with a value at the Moon’s centre of 
500° K. 


where B — 0.9 X 10“^ cal cm-^ s-^ considers 
heating the Moon by solar radiation. On cal- 
culation this condition leads to practically the 
same results as does the more simple condi- 
tion T = 273° K = const. 

THE MECHANISM OF DIFFERENTIATION 

The mechanism of differentiation is de- 
fined by the function for the distribution of 
radioactive elements F{r,t) which is given 
as follows. At first, we considered a homo- 
geneous model for the distribution of radio- 
active elements, i.e., F (r,t) = 1 (analogous 
to refs. 22 and 23) for the first stage of the 
calculation. Then, after reaching at the 
moment U, at the point n, the temperature 
corresponding to the beginning of melting 
Tinit, melt, the distribution of temperature ob- 
tained up to the time t 2 = fi -H At (where 
At is the time of melting (see below) is 
taken as the new initial distribution, and the 
distribution of radioactive elements is as- 
sumed to be layered ; namely, radioactive ele- 
ments from the layer ± Ar are transferred 
to the adjacent upper layer (r£ ± Ar) and 
the calculation is carried out with the de- 
pleted layer (ri ± Ar). However, for the 
time At defined by the necessity of heating 
the layer by AT = 200° K, deeper layers have 
reached the melting temperature (since 
Tinit. melt, increases with the depth, i.e., with 
the decrease of r) . Radioactive elements 
from these layers go to the layer ri zt Ar, i.e., 
their content in this layer is equalized and so 
on, until the lower layers no longer melt. 
Such layered differentiation (similar to that 
used by Fricker, Reynolds, and Summers) is 
defined by five variants that differ in the 
value n which is the fraction of radioactive 
elements removed from the layer. We accept 
n — \ (complete removal), 0.8 (80 percent 
of the radioactive elements are removed from 
the layer), 0.6, 0.4, and 0.2. In this way, 
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Figure 1 . — The time dependence of temperature in different layers of the Moon for n — 1 (along the radius, 
starting from the center, 0, layer SO — S56 km, i0,60 — 96 — 1709 km). The melting temperature corres- 
ponding to each layer is shovm (dashed lines), as is the time dependence of the start of melting (heavy 
solid line). The times of initial melting of the Moon, the end of melting of the outer layers, as well as 
the start of differentiation are shown on the abscissa. 


after differentiation in each of the melted 
layers, there remains from 0 to 80 percent 
of the content of radioactive elements in the 
non differentiated material (note that in the 
latest paper by Mayeva (ref. 13) the residue 
was from 2 to 20 percent, i.e., n >0.8). 

The present paper does not consider the 
possibility of the melted material’s penetrat- 
ing into the external solid layer as was done 
in the paper by Mayeva (ref. 13), since the 
concentration of radioactive elements in the 
thin surficial layer leads, in reality, to more 
rapid cooling of the Moon and to the greater 
thickness of the hard crust, as obtained in 
the model of exponential differentiation 
considered in detail in references 22 and 23. 

The Methofd of Calculation 

For the numerical solution of the heat 
conductivity equation we use the method of 
line coordinates where the section [0,ro], 
(where „ is the radius of the Moon) is divided 


into 100 layers and only the differentials in 
the spatial coordinate were replaced by dif- 
ference relations ratio. Thus the solution of 
the boundary problem is reduced to a solution 
of the Cauchy problem for the system of ordi- 
nary differential equations. It should be noted 
that the approximation error is larger, the 
closer the layer is to the surface. That is why 
a non-uniform grid is used in the calculation : 
the layer thickness close to the surface («* 5 
km) is several times smaller in comparison 
with the depth layer (*=« 18 km). The system 
obtained is solved by the approximation 
method, which provides an accuracy of the 
temperature calculation (at each step in 
time) of the order of 10~®. All calculations 
are made by the computer BESM-6. 

Results of Calculation 

Results of the calculation are given in fig- 
ures 1 through 12 and in tables 2 and 3. 
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Figure 2 . — The time dependence of temperature in different layers of the Moon. 


THE PROCESS OF MELTING AND 
DIFFERENTIATION 

Figures 1 through 4 show the time depen- 
dence of temperature in different layers of 
the Moon (along the radius beginning from 
the centre 0) for all variants at n = 1. Tem- 
peratures at the beginning of melting are 
plotted for each layer and the time depen- 
dence of the beginning of melting in different 
layers. From figures 1 through 4 it is seen 
that melting begins in layer 80^ for all 
variants, about 1500 km from the center (i.e., 
at the depth 250 km under the surface) 


'Altogether there are 100 layers. 


during the period from 0.4 to 1.5 billion years 
from the beginning of lunar history. Rapidly 
(in Oil to 0.3 billion years) the region of 
melting propagates to the center and to the 
Moon’s surface. In the layer where melting 
first begins, differentiation of the material 
begins after 0.2 to 0.5 billion years. The re- 
gion of melting at the start of differentiation 
(tinit. dif.) reaches 1600 km; at the same time 
there is a sharp temperature increase in the 
layers close to the surface. As a result the 
temperature in these layers is close to or ex- 
ceeds Tinit. melt.. This is Understandable since 
the melted radioactive elements are moved 
by convection to these layers from deeper 
layers. The deep layers (at a depth of more 
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than 500-600 km, i.e., r = 1100-1200 km) 
prior to differentiation have been heated 
above T^t. meit. and are now slowly cooling 
due to the convection of radioactive elements 
to the upper layers^. From figures 5 through 
7 it is seen that in this way almost simul- 
taneous heating of lunar material occurs up 
to depths of 200-250 km under the surface. 
Layers which are closer to the surface heat 
more slowly and the thickness of these layers 
increases with time. Numerical data are given 
in tables 2 and 3. 

“ Differentiation occurs in all variants except for 
III (cold model) , where the energy of radioactive 
elements is not sufficient for the beginning of 
differentiation. 


THE BEGINNING AND THE DURATION 
OF THE MAXIMUM MELTING OF THE 
INTERIOR 

Figures 8 through 11 give the time depen- 
dence for the radius of melting. From these 
figures it is seen that in all variants the pe- 
riod of the maximum melting continues for 1 
to 2 billion years ( depending on the variant 
— see tables 2 and 3) and begins in the inter- 
val from 0.7 to 2.3 billion years from the 
beginning of lunar history. Naturally, for 
large concentrations and “hot” models the 
maximum melting begins earlier than for low 
concentrations and “cold” models. So for the 
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Figure 4 . — The time dependence of temperature in different layers of the Moon. 


variants I and the maximum melting be- 
gins approximately 3.5 to 4 billion years 
prior to the present time and continues for 
1.5 to 1 billion years. In the case of low 
(variant III) and intermediate concentra- 
tions and large convection (n > 0.8) of ra- 
dioactive elements (variant II) the maximum 
melting begins from 2.5 to 3 billion years 
prior to the present and continues for 1.5 to 
2 billion years, i.e., it begins later and con- 
tinues longer, as expected*). A special case is 
variant II (“cold” model), which has a low 
degree of convection of radioactive elements 

* The smaller the concentration, the more quietly 
the processes occurred. 


(n < 0.6). In this case the period of max- 
imum melting begins 2.5 billion years prior 
to the present and continues to the pres- 
ent time. Here a specific combination of 
conditions occurs. The concentration of ra- 
dioactive elements is not so large and the 
process of differentiation is not so effective 
that a fast convection of radioactive elements 
to the surface occurs. Thus their influence 
on further history is small, but not suffici- 
ently small for the natural cooling of the 
Moon. 

Consequently, in all variants except for 
the above-mentioned II (n < 0.6) the Moon 
is presently cooling. 
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Figure 5 . — The temporal variation of temperature in the lunar interior. 


THE THICKNESS OF THE CRUST 

As it is seen from figures 8 through 11 and 
tables 2 and 3, the thickness of the crust at 
the period of maximum melting ranges in 
variants I, II, and III from 15 to 20 km, 
from 25 to 30 km, from 35 to 45 km, and 
from 100 to 200 km, respectively; i.e., from 
1 to 2.5 billion years ago the Moon was almost 
completely melted (in variant II with n < 
0.6, the Moon even now should be melted up 
to the depth of about 45 km under the sur- 
face). Then cooling began and presently the 
crustal thickness in all variants is from 150 


to 250 km (except for the above-mentioned 
“special case”) . The same thing follows from 
figure 12, where the temperature distribution 
is given in the Moon for the present time. 


THE HEAT FLUX THROUGH THE LUNAR 
SURFACE 

Figure 12 gives the time dependence of the 
heat fiux through the lunar surface. From 
these curves it is seen that at the period of 
maximum melting of the Moon the flux, as 
could be expected, is also maximum and ex- 
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Figure 6 . — The temporal variation of temperature in the lunar interior. 


ceeds 2 to 3 times the stationary value for 
the given concentration of radioactive ele- 
ments. At the present time for all variants 
(except for II with n < 0.6), the flux is de- 
creasing and is approaching the stationary 
value 9stat., yet even now the flux is nearly 1.5 
times 9atat.- The flux value in variants I, I+% 
II {n > 0.8), and III is equal to 1.8 X 10"® 
cal cm'2 S"S 1.27 X 10 ® cal cm VS and 0.9 
to 0.95 X 10"® cal cm"* s"^ respectively. In 
variant III (“cold” and “hot” models) it is 
0.28 and 0.55 X 10 ® cal cm"* s"^ and variant 
II (n < 0.6) has a rather large flux equal to 
1.6 X 10 ® cal cm * s"^, which is associated 


with the continuation of melting in this vari- 
ant up to the depth of 45 km. 

Thus fluxes close to those measured by 
radio astronomical methods and obtained by 
Apollo 15 are realized in variants and II 
(n > 0.8), i.e., for the terrestrial or “man- 
tle” pyrolite compositions of lunar material. 

Conclusion 

Calculations are given for different con- 
centrations of radioactive elements corres- 
ponding to a chondritic composition (III), to 
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a terrestrial composition (II), and a terres- 
trial mantle-pyrolite composition (I and I^). 
Differentiation of material and the convec- 
tion of radioactive elements to the surface 
due to melting are taken into account in the 
calculation by a scheme similar to that used 
by Fricker, Reynolds, and Summers. The 
calculations show that the Moon’s interior 
has been heated to melting during the first 
0.7 to 2.3 billion years after its formation. 
The maximum melting of the Moon covered 


practically the entire Moon up to a depth of 
15 to 45 km below its surface and began from 
3.5 to 4 billion years ago, if one assumes a 
pyrolitic composition of the Moon; or 2.5 to 
3 billion years ago, if one assumes a terres- 
trial or chondritic composition of the Moon. 
This maximum melting continued for 1 to 2 
billion years and at present the Moon is 
cooling (true, there is an improbable variant 
of weak differentiation for the “terrestrial” 
composition of the Moon when melting is 
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occurring even now). The present thickness 
of the crust is 150 to 250 km and the heat 
flux exceeds the approximately 1.5 times sta- 
tionary value. The flux value for variants I, 
II, and III is equal to 1.8, 1.27, 0.9 to 
0.95, and 0.28 to 0.55 X 10-* cal cm'^ s'*, re- 
spectively. 

From tables 2 and 3 it follows that in the 
most probable variants with the pyrolite- 
“mantle” (P+) and terrestrial (II) composi- 
tion of the initially “hot” Moon, the 


maximum melting of the Moon began 3.3 to 
3.8 billion years ago, which corresponds to 
the age of mare basalts and mascons (3.16 
to 3.71 billion years) . These data testify to the 
volcanic activity that was completed aproxi- 
mately 2.7 to 3 billion years ago — evidence in 
agreement with the absence of lunar samples 
younger than 3 billion years (refs. 32, 33, 
and 34), when the activity stopped. Thus, 
limitations stated recently by Toksoz and 
Solomon (ref. 21) are satisfied, yet the thick- 


362 


COSMOCHEMISTRY OF THE MOON AND PLANETS 



Figrure 11. — The distribution of present temperatures in the Moon as a function of n. 


ness of the crust obtained by our calculations 
is essentially small (ref. 21) — a crust 600 km 
thick. 

It should be noted that the most realistic 
variant uses a “terrestrial” composition for 
the Moon (composition “C” according to 
Starkova), which gives, in spite of its mod- 
erate concentration of radioactive elements, 
a sufficiently large heat flux^ that corres- 
ponds to the radio astronomical measure- 
ments and Apollo 15 data, i.e., to obtain the 
large observed flux it is not necessary to have 
high concentrations of radioactive elements. 

' This possibility was noted by Lubimova in ref- 
erence 6. 
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Table 3 . — Calculated Results for Various Thermal Model Variants 



(cal cm'* sec'*) lO"* 


tm.i mel flux 

rmelt init dif*” 

Ar crust min 

Ar” 
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Note: (1) n.it m.it and ri.it utt = 1500 km. 
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